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Abstract 

Expressing explicitly the parameters of the standard Skyrme interaction in terms of the macro- 
scopic properties of asymmetric nuclear matter, we show in the Skyrme-Hartree-Fock approach 
that unambiguous correlations exist between observables of finite nuclei and nuclear matter prop- 
erties. We find that existing data on neutron skin thickness Ar^p of Sn isotopes give an important 
constraint on the symmetry energy Egymipo) and its density slope L at saturation density po- 
Combining these constraints with those from recent analyses of isospin diffusion and double neu- 
tron /proton ratio in heavy- ion collisions at intermediate energies leads to a more stringent limit on 
L approximately independent of -EsymCPo)- The implication of these new constraints on the Ar^p 
of ^"^^Pb as well as the core-crust transition density and pressure in neutron stars is discussed. 

PACS numbers: 21.65.Ef, 21.10.Gv, 26.60.Gj, 21.30.Fe 
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I. INTRODUCTION 



The nuclear symmetry energy Esym{p) that encodes the energy related to the neutron- 
proton asymmetry in the equation of state (EOS) of isospin asymmetric nuclear matter 
(ANM) plays a crucial role in both nuclear physics and astrophysics jl-6|. It is also relevant 
to some interesting issues regarding possible new physics beyond the standard model j^- 



10|. Although significant progress has been made in recent years in determining the density 



dependence of -E'sym(p) p, ISj, large uncertainties still exist even around the normal density 
pQ. For instance, while the value of -EsymlPo) is determined to be around 30 ± 4 MeV mostly 
from analyzing nuclear masses, the extracted density slope L of Esymip) at po scatters in 
a very large range from about 20 to 115 MeV depending on the observables and methods 
used in the studies lll-ll3l|. Since many observables in terrestrial laboratory experiments 



intrinsically depend on both -Esym(Po) and L, the extraction of L at an accurac y r equired 
for understanding more precisely many important properties of neutron stars jsl-ll^ is still 
severely prohibited, although the uncertainty of -E'sym(Po) is relatively small. To extract L 
with higher accuracy is thus of crucial importance. 

heoretically, studies based on both mean-field theories 15N21| and droplet-type models 

^p)^/^ of heavy nuclei, 



22 



25| have shown that the neutron skin thickness Ar, 



np 



2\l/2 



given by the difference of their neutron and proton root-mean-squared radii, provides a 
good probe of Esy^iip)- In particular, Ar„p has been found to correlate strongly with both 
-Esym(po) and L in microscopic mean-field calculations |l5l-l2l| using different parameter sets 
for the nuclear effective interactions, which all fit the binding energies and charge radii of 
finite nuclei but correspond to different £'sym(p) and give different Ar„p. It is, however, 
difficult to extract an accurate value for L from comparing calculated Ar„p of heavy nuclei 
with experimental data as Ar„,p depends on a number of nuclear interaction parameters in 



a highly correlated manner 



17, 



18l ] and the calculations have been usually carried out by 



varying simultaneously the interaction parameters. Similar difficulties also exist when one 
tries to extract other physical quantities from observables of finite nuclei within mean-field 



theories or density functional theories 26|]. A well-known example is the Skyrme-Hartree- 



Fock (SHF) approach using normally 9 interaction parameters. Although experimental data 



on nucleon-nucleon scatterings and properties of both finite nuclei and infinite nuc 
would in principle put strong constraints on the combinations of these parameters 



ear matter 



19 



27 



32|, 
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there is generally no constraint on most of the individual interaction parameters. Instead of 
varying directly the 9 interaction parameters within the SHF, we propose here an alternative 



approach based on a modified Skyrme-Like (MSL) model [33|] to express them explicitly 
in terms of 9 macroscopic observables that are either experimentally well constrained or 
empirically well known. This opens the possibility to explore transparently the correlations 
between properties of finite nuclei and the macroscopic properties of nuclear matter within 
the SHF approach. In the present work, we use this method to study the correlation between 
Ar„p and various macroscopic observables of infinite nuclear matter by varying individually 
the values of the latter within their known ranges. We then demonstrate that existing Ar^p 
data on Sn isotopes can give important constraints on L and -E'sym(Po)- Combining these 
constraints with those from recent analyses of isospin diffusion and double neutron/proton 



ratio in heavy-ion collisions at intermediate energies 



we further show that a more 



stringent limit on L is obtained approximately independent of the value of EsymiPo)- Finally, 
we discuss the implication of these new constraints on both the Ar„p of ^°'^Pb and the core- 
crust transition density and pressure in neutron stars. 



II. THE THEORETICAL MODEL 



In the present work, we use the so-called standard Skyrme interaction (see, e.g., Ref. 34|). 



which has been shown to be very successful in describing the structure of finite nuclei, 
especially global properties such as binding energies and charge radii, although non-standard 
extension is possible js^. In the standard SHF model, the total energy density of a nucleus 



is written as 



3^ 



v. — IC + T-Lq + T-Ls + T-Leff + 'Hfin + 'Hso + ^sg + 'Hcoul (1) 
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where /C = is the kinetic-energy term and T-Lcoui is the Coulomb term, and "Ho, "Hs, 
'Heff, Ufin, Hso, 'Hsg are given by 

-Ho = to[(2 + Xo)p' - (2x0 + 1){pI + pI)]/^ (2) 
= hp''[{2 + xs)p' - (2x3 + 1)(p5 + p')]/24 (3) 
Heff = [t2(2x2 + 1) - ti(2xi + l)](r„p„ + rpPp)/8 

+ [ti(2 + Xi)+t2(2 + X2)]rp/8 (4) 
Ufin = [3ti(2 + xi) -t2(2 + X2)](Vp)V32 
-[3ti(2xi + l)+t2(2x2 + l)] 

X [(Vp„)^ + (Vpp)^] /32 (5) 

= Wolf ■^P + ^p-^Pp + ^n-^Pn]/2 (6) 

Hsg = (ti-t2)[Jj + J^]/16 

-(tlXi +t2a;2)JVl6 (7) 

in terms of the 9 Skyrme interaction parameters a, to ~ ^3? — ^3, and the spin-orbit 
couphng constant Wq. In the above equations, pi, Ti and Jj are, respectively, the local 
nucleon number, kinetic energy and spin densities, whereas p, r and J are corresponding 
total densities. 

In the MSL model, the EOS of symmetric nuclear matter (SNM) and the nuclear sym- 
metry energy i?sym(p) can be expressed, respectively, as [J 



Eo{p) = E^,y/^ + + au/2 + /3uV(7 + 1), (8) 

EsyM = E^Z.{poW" + Du'/' + E^Up) (9) 



where u = p/ po is the reduced density; E^^^ and E^^ are, respectively, the kinetic energy 
at po and its contribution to £'sym(p); and El°'^{p) is the local density-dependent symmetry 
energy given by 

E^M = (1 - y)E':;Mu + yEl;i,{poW (lo) 

with the dimensionless parameter 

L - 3E,,M + E't^M - 2D 
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The model also includes following density-gradient term in the interaction part of the binding 
energies for finite nuclei 



^grad = Gsiypf/{2p) - Gv [V(p„ - p,)f /(2p), 



(12) 



where Gs and Gy are the gradient and symmetry-gradient coefficients 28 1. 

By comparing expressions ([8]), (|9]), and f|T2|) in the MSL model with corresponding ones in 
SHF, the 9 Skyrme interaction parameters in Eqs. ©-([Tj) can be related to the 9 parameters 
a, /9, 7, C, £'g°^(po), 2/, Gs and Gy in the MSL model by following analytic relations 



to 

Xq 

ts 

tl 
t2 

Xi 

a 



4«/(3po) 

3(y- 1)^1- (po)/«- 1/2 

16/3/ (7 + 1)] 

-3y{^ + l)El;M/{2(5)-l/2 

20C/ [9poi4y] + 8G5/3 
4(25C - 18D) 8{Gs + 2Gv) 



9po(fc^)2 



12Gv 



6D 
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20Gv + 4G5 - 
1 



Po(4) J 
5(16C- 18D) 

Spoik'^y 



/(3t2 



(13) 
(14) 
(15) 
(16) 
(17) 
(18) 

(19) 

(20) 
(21) 



with kp 



I.Stt^Po) • Since the 7 parameters a, /3, 7, C, D, _E''°^(po) and y in the MSL 
model can be expressed analytically in terms of the 7 macroscopic quantities po, -E'o(po), 
ity Kq, the isoscalar effective mass m*Q, the isovector effective mass m*g, 



33|, the 9 Skyrme interaction parameters a, to — ^3, Xq — X3 can also be 



the incompressibi 

^sym(Po), and L 

expressed analytically in terms of the 9 macroscopic quantities po, -Eo(po), Kq, m*^, m* Q, 
EsyraiPo): L, Gs, and Gy via above relations. 

As a reference for the correlation analyses below, we use following default values for 
the macroscopic quantities. For the 7 bulk properties of ANM, we take po = 0.16 fm~^, 
Eo(po) = -16 MeV, Kq = 230 MeV, m*Q = 0.8m, m^^ = 0.7m, Ssym(po) = 30 MeV, 
and L = 60 MeV. Empirically, values of the gradient coefficient Gs and the symmetry- 
gradient coefficient Gy are poorly known, although they can be constrained by the nuclear 
surface energy coefficient as, the Landau parameters for the spin and spin-isospin channels 
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FIG. 1: (Color online) The symmetry-gradient coefficient Gy dependence of Gs and VFq (a), 
Landau parameters Go(po)5 G^oCpo)) G'i(po)) Go(po) (b), and the reduced critical density pcr/ Po for 
symmetric nuclear matter and pure neutron matter (c). 

in symmetric nuclear matter at the saturation densit y, i .e., (jn(Pn), CqIpo)) CiIpo) and 



31 



(j"q(po); and the stabihty condition of nuclear matter 
Skyrme parameter sets, we have roughly Gs = 110 ~ 
MeV-fm^. Here we use the empirical value of the surface energy coefficient as = 18 MeV 22 1 



32] . For the existing standard 
150 MeV-fm^ and Gy = -40 ~ 40 



to determine Gs- We note t 
orbit coupling constant Wq 



lat the surface energy coefficient as also depends on the spin- 



31 



32 



35| that can be determined by the neutron pi/2 — P3/2 
splitting in ^^O. Since the latter depends on Gs and Gy, the three quantities Gs, Gy, and 
Wq need to be determined simultaneously. 

Keeping the 7 bulk properties of ANM unchanged, we plot in Fig. [1] (a) Gs and Wq as 
functions of Gy with the values of Gs, Wq, and Gy simultaneously giving the surface energy 
coefficient as = 18 MeV and fitting the neutron pi/2 —P3/2 splitting in ^^O. In Figs. [1] (b) and 
(c), we further show the Gy dependence of the Landau parameters Go{po), Gq{po), Gi(po), 
G'q{pq), and the critical density pcr above which at least one Landau parameter violates the 
stability condition for symmetric nuclear matter and pure neutron matter. Requiring the pcr 
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TABLE I: Skyrme parameters in MSLO (left side) and some corresponding nuclear properties (right 
side). 



Quantity 


MSLO 


Quantity 


MSLO 


to (MeV-fm^) 


-2118.06 


On ffm""') 


0.16 


h CMeV-fm^l 


395 196 

KJ %J \J • i \J 


En (MeV) 


-16.0 


t2 (MeV-fm^) 


-63.9531 


Kn (MeY) 


230.0 


is (MeV-fm^+^'^) 


12857.7 




0.80 


Xq 


-0.0709496 


ml n/m 


0.70 


Xl 


-0.332282 


^sym(po) (MeV) 


30.0 


X2 


1.35830 


L (MeV) 


60.0 


X3 


-0.228181 


Gs (MeV-fm^) 


132.0 


a 


0.235879 


Gv (MeV-fm^) 


5.0 


Wo (MeV-fm^) 


133.3 




0.42 



of symmetric nuclear matter to be larger than po leads to Gy ^ —20 MeV-fm^ and putting 
Per > Po for pure neutron matter further leads to Gy ^ 5 MeV-fm^. Empirically, the Landau 
parameter G'q has been extensively investigated, and its value can vary from about zero to 



1.6 depending on the models and methods 



32 



36 



- |42| . In the present work with the standard 



SHF approach, a positive G'q leads to Gy ^ 70 MeV-fm^ as shown in Fig. [T](b). Furthermore, 
the Landau parameter G'q can be extracted from the spin-isospin response in finite nuclei, 
and its value has been found to be 0.45 ± 0.06 in the standard SHF approach 43|, |44 |. 
Therefore, we choose here Gy = 5 MeV-fm^ which leads to G'q = 0.42, Gs = 132 MeV-fm^, 



ue Gs = 132 MeV-fm^ is 



45( 1 . This new Skyrme 



and Wq = 133.3 MeV-fm^. It is interesting to see that the va 
quite consistent with that used extensively in the literature fl^ . 
parameter set obtained with above empirical values for the macroscopic quantities is referred 
as MSLO. Summarized in Table [T] are values of corresponding Skyrme parameters and some 
macroscopic quantities. 

To test the new Skyrme parameter set MSLO, we calculate the binding energies and 
charge rms radii for a number of closed-shell or semi-closed-shell nuclei: ^^O, ^°Ca, ^^Ca, 
^^Ni, ''^Ni, ^°Zr, ^°°Sn, ^^^Sn, and ^°®Pb. Figure [2] shows the relative deviation of the charge 



rms radii and binding energies of these nuclei from those measured in experiments 46N48|. 
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FIG. 2: (Color online) Relative deviation of the binding energies and charge rms radii of ^^O, ^'^Ca, 
48Ca, s^Ni, 78Ni, 90Zr, ^oosn, ^^^Sn, and ^ospb from SHF with MSLO. 

It is seen that the MSLO can describe the experimental data very well except for the light 
nucleus ^^O for which the deviation reaches to about 2 — 3%. This is a remarkable result 
as MSLO is not obtained from fitting measured binding energies and charge rms radii of 
finite nuclei as in usual Skyrme parametrization. It should be pointed out that our main 
motivation for introducing the MSLO is not to construct another Skyrme parameter set to 
describe data, but to use as a reference for the correlation analyses in the following. As we 
will show, varying Gs, Gy and will not affect the conclusion in the present work. 



III. RESULTS 



To reveal clearly the dependence of Ar„p on each macroscopic quantity, we vary one 
quantity at a time while keeping all others at their default values in MSLO. Shown in Fig. 
|3]are the values of Ar„p for ^°^Pb, ^^°Sn and ^^Ca. Within the uncertain ranges considered 
here, the Ar„p of ^^^Pb and ^^°Sn exhibits a very strong correlation with L. However, it 
depends only moderately on Esym{po) and weakly on m*Q. On the other hand, the Ar„p 
of ^^Ca displays a much weaker dependence on both L and -Esym(po)- Instead, it depends 
moderately on Gy and Wq. This explains the weaker Ar„p-i?sym(p) correlation observed for 
^^Ca in previous SHF calculations using different interaction parameters 20|. These results 
demonstrate that the Ar„p of heavy nuclei can provide reliable information on the symmetry 
energy around the normal density. 

As we vary one of the macroscopic quantities in Fig. |3l it is of interest to see how this 
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FIG. 3: (Color online) The neutron skin thickness Ar„p of ^"^^Pb, ^^"^Sn and ^^Ca from SHF with 
MSLO by varying individually L (a), Gy (b), Gs (c), ^o(po) (d), £^sym(/Oo) (e), (f), m*o (g), 
"^^,0 (h), Po (i), and (j). 

affects the binding energy and charge rms radius. This is shown in Fig. H] and Fig. [5l 
respectively, for the changes in the relative deviation of the binding energies and charge 
rms radii of ^°^Pb, ^^°Sn and "^^Ca from the data. Within the uncertain ranges for the 
macroscopic quantities considered here, it is seen that the relative deviation of the binding 
energy is basically less than 3% except the case for the macroscopic quantity E^. For the 
charge rms radius, the relative deviation is even much smaller. Especially, for the heavy 
^°^Pb and ^^"^Sn, the relative deviation of the charge rms radius is basically less than 0.5% 
except the case for the macroscopic quantity po- These features imply that the binding 
energies and charge rms radii of finite nuclei can be reasonably reproduced when we perform 
the correlation analysis shown in Fig. [3] by varying individually the macroscopic quantity. 
This is particularly the case when we only change the macroscopic quantities L and -E'sym(Po) 
from the MSLO in the correlation analyses. 

Experimentally, much effort has been devoted to determining the values of Ar„p for finite 
nuclei using various methods. In particular, the Ar„p of heavy Sn isotopes has been sys- 
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FIG. 4: (Color online) Relative deviation of the binding energies of ^'^^Pb, ^^'^Sn and ^^Ca from 
SHF with MSLO by varying individually L (a), Gy (b), Gs (c), So(po) (d), ^sym(po) (e), (f), 
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FIG. 5: (Color online) Same as Fig. U] but for the charge rms radii. 
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FIG. 6: (Color online) (a): The Ar„p data for Sn isotopes from different experimental methods 
and results from SHF calculation using MSLO with L = 20, 60 and 100 MeV. (b): & function 

of L. (c): Constraints on L and Esymipo) from the analysis of the Ar„p data on Sn isotopes 
(Grey band as well as dashed and dotted lines). 



tematically measured 49N54j. As an illustration, we first show in Fig. [6] (a) the comparison 
of available Sn Ar„p data with our calculated results using different values of 20, 60 and 
100 MeV for L and the default values for all other quantities in MSLO. It is seen that the 
value L = 60 MeV best describes the data. To be more precise, the evaluated from the 
difference between the theoretical and experimental Ar„p values is shown as a function of 
L in Fig. [6] (b). The most reliable value of L is found to be L = 54 ± 13 MeV within a 2a 
uncertainty. 

Since the value of Ar„p depends on both L and -E'sym(Po)! ^ two-dimensional analysis 
as shown by the grey band in Fig. [6] (c) is necessary. It is seen that increasing the value 
of -E'sym(po) systematically leads to smaller values of L. More quantitatively, the value of 
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FIG. 7: (Color online) Contour curves in the Esym{po)-L plane for the Ar„p of ^"^Pb (a) from SHF 
calculation with MSLO and the core-crust transition density pt (b). The shaded region represents 
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the overlap of constraints obtained in the present work (dashed lines) and that from Ref. 
(dash-dotted lines). 

L varies from 67 ± 10.5 to 37 ± 15.5 MeV when the value of -E'sym( Po) changes from 26 
to 34 MeV. Furthermore, we have estimated the effects of nucleon effective mass by using 



m*Q = 0.7m and m*Q = 0.6m as wef 



the empirical constraint m*Q > m*Q j6,[55|, and the resulting constraints are shown by the 
dashed and dotted lines. As expected from the results shown in Fig. |3l effects of nucleon 
effective mass are small with the value of L shifting by only a few MeV for a given i?sym(po)- 
We have also checked that effects of varying other macroscopic quantities are even smaller. 

The above constraints on the L-i?sym(po) correlation can be combined with those from 
recent analyses of isospin diffusion and double n/p ratio in heavy- ion collisions at interme- 
diate energies [111 to determine simultaneously the values of both L and i?sym(po)- Shown 
in Fig. [7] (a) (and (b)) are the two constraints in the -Esym(po)--^ plane. Interestingly, these 
two constraints display opposite L-£'sym(Po) correlations. This allows us to extract a value 
of L = 58 ± 18 MeV approximately independent of the value of -E'sym(po)- This value of L is 
essentially overlapped with other constraints extracted from different experimental data in 
131] but with much higher precision although the constraint on -E'sym(Po) 



as m*Q = 0.9m and m*Q = 0.8m, in accord with 



the literature 
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is not improved. It also agrees well with the value of L = 66.5 MeV obtained from a recent 
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FIG. 8: Same as Fig. [3] but for the core-crust transitfon density pt in neutron stars. 



systematic analysis of the density dependence of nuclear symmetry energy within the micro- 
scopic Brueckner-Hartree-Fock approach using the realistic Argonne V18 nucleon-nucleon 
potential plus a phenomenological three-body force of Urbana type 56|. Furthermore, it 
is in remarkably good agreement with the value of L = 52.7 MeV extracted most recently 
from global nucleon optical potentials constrained by world data on nucleon-nucleus and 
(p,n) charge-exchange reactions 57 1. 

Also shown in Fig. [7] (a) are contours of the Ar„,p of ^"^^Pb. Based on the constraints 
on L and -E'sym(po) shown by the shaded region in Fig. [3, it is seen that the Ar„p of ^osp]-, 
is tightly limited to a narrow region of 0.175 ± 0.02 fm, which is quite consistent with 
other constraints from various 

Lead Radius Experiment (PREX) 58|] being preformed at Jefferson Lab aims to determine 
model- independently the (r^)^/^ of ^osp]^ i^y^ accuracy, and this is expected to further 
improve the determination of Esy^{p) at subnormal densities. 

To see the implications of our results in astrophysics, we have carried out a similar 
correlation analysis for the transition density pt and corresponding pressure Pt at the inner 
edge of neutron star crusts, which play crucial roles in neutron star properties , usmg 

their values evaluated in a dynamical approach [l^ , and the results are shown in Fig. [8] and 
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FIG. 9: Same as Fig. [3] but for the the core-crust transition pressure Pt in neutron stars. 



Fig. [9l respectively. It is seen that the pt (Pt) displays a particularly strong correlation with 
L (L and Egymipo)), a weak dependence on i?sym(po) and Kq (Kq), but almost no sensitivity 
to other macroscopic parameters. These features are consistent with the results in Ref. 14 1 
where pt has been shown to display a stronger correlation with L than Pt in SHF calculations 
using different interaction parameters. The contours of the core-crust transition density pt 
in neutron stars in the Esym{Po)-L plane is shown in Fig. [7| (b). It shows that the value of 
Pt is limited to 0.069 ± 0.011 fm~^ by the constraints on L and £'sym(po) obtained in the 
present work. Including further the uncertainty in the value of Kq, we obtain a value of 
Pt = 0.069 ± 0.018 fm-3. A similar analysis leads to Pt = 0.33 ± 0.21 MeV/fm^. These 
results agree well with the empirical values but are slightly larger than previous results 



in Ref. 



141 ] using Esymipo) = 30.5 MeV and L = 86 ± 25 MeV extracted only from the 



isospin diffusion data in heavy- ion collisions [59|. 



IV. SUMMARY 



We have proposed to analyze the correlation between observables of finite nuclei and some 
macroscopic properties of asymmetric nuclear matter by expressing explicitly the parameters 
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of the nuclear effective interaction in terms of the macroscopic properties of asymmetric 
nuclear matter. This would allow us to extract information on some important physical 
quantities from data on finite nuclei in a more transparent way. 

Using such a correlation analysis within the standard SHF approach, we have demon- 
strated that the neutron skin thickness of heavy nuclei can provide reliable information on 
the symmetry energy, and the existing neutron skin data on Sn isotopes can give important 
constraints on the symmetry energy parameters Esymipo) and L. In particular, combining 
the obtained L-Esymipo) constraints with that from recent analyses of isospin diffusion and 
double n/p ratio in heavy-ion collisions has led to a quite accurate value of L = 58± 18 MeV 
approximately independent of the value of -E'sym(Po)- The obtained L-EsymiPo) constraints 
also put a stringent limit of Ar„p = 0.175 ± 0.02 fm for the neutron skin thickness of ^"^^Pb. 

Furthermore, we have explored how the core-crust transition density pt and the corre- 
sponding pressure Pt in neutron stars correlate with the macroscopic properties of asym- 
metric nuclear matter. Our results have indicated that the pt displays a particularly strong 
correlation with L, a weak dependence on Esym{po) and Kq, but almost no sensitivity to 
other macroscopic parameters. On the other hand, the Pt exhibits a strong correlation with 
both L and Esymipo), and a weak dependence on Kq. The L-E^y^^po) constraints obtained 
in the present work leads to pt = 0.069 ± 0.018 fm^^ and Pt = 0.33 ± 0.21 MeV/fm^. 

Although we have mainly analyzed in the present work the correlation between the neu- 
tron skin thickness of finite nuclei and the symmetry energy within the standard SHF ap- 
proach, our method can be generalized to other correlation analyses or mean-field models. 
In particular, it will be interesting to see how our results will change if different energy 
density functions are used. These studies are in progress. 
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